Low-energy potential scattering in two and three dimensions 



O 

o 

(N 

c 
3 



N. N. Khuri/'Q Andre Martin, 2 ^ J.-M. Richard, 3 '^ and Tai Tsun Wu 4 ' 2 f] 

Department of Physics, The Rockefeller University, New York, New York 10021 
2 TH Division, CERN, CH - 1211 Geneva 23 
3 Laboratoire de Physique Subatomique et Cosmologie, 
IN2P3-CNRS, Universite Joseph Fourier, INPG, 

53, avenue des Martyrs, Grenoble, France 
^ Gordon McKay Laboratory, Harvard University, 
Cambridge, Massachusetts 02138-2901 
(Dated: June 10, 2009) 

Conditions are established for the existence of a scattering length and an effective 
range in the low-energy expansion of the S-wave phase-shift of a central potential 
in two and three dimensions. The behavior of the phase-shift as a function of the 
momentum is also derived for longer-range power-law potentials which do not fulfill 
these conditions. 
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INTRODUCTION 



Ten years ago, low-energy scattering in two spatial dimensions (2D) was studied for a 
large class of potentials [lj. It was found that the S-wave phase-shift as a function of the 
momentum k, 5o(k), has a universal behavior 



tan 5 



7[ 



1 



0(1) 



2 In k Ink 

as k — > 0. See, e.g., [2JE1I1] for earlier studies, and [3] for the case of non-central potentials. 

This behavior is radically different from that for three spatial dimensions (3D), where 
5o(k) ~ — ak, where a is the scattering length. 



This universal behavior (1.1) is valid for most potentials in 2D, but there is a subclass of 



potentials where Eq. (1.1) does not hold. Special cases of (1.1) have been studied earlier [T]. 

This universal behavior is the underlying reason why wire antennas are very efficient and 
widely used [6]. 



It is difficult to obtain an approximate value for tan 5 from this asymptotic formula (1.1). 
The basic reason is that, since k has the dimension of inverse length, the value of In k depends 
on the length unit used. In order to get a good approximate value, we need some knowledge 
of the next to leading term, i.e., we would like to rewrite the above result in the form 

1 



tan 5 



71 
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for small k, and obtain more information about the scale ho- What can ko be? Clearly 
it cannot be universal. And it is reasonable to think that the value of k is related to a 
scattering length. 

We recall that in 3D, the scattering length appears in two different ways. For small k, we 
have 8 ~ — ka (with the convention that a > for a repulsive potential). But in addition, 
the scattering length, a, can be defined by the zero of the reduced zero-energy asymptotic 
wave- function 

Masym = C ' (r - a) , (1.3) 

We use this latter route to define a scattering length in 2D. We shall see that the scale in 



(1.2) is proportional to the inverse of a, the value for which the zero-energy asymptotic wave 



function vanishes. 



In Sec. II, we treat, for pedagogical reasons, the simple case of a potential with finite 



range. In Sec. Ill, we derive sufficient conditions on the potential for the existence of the 



scattering length in 2D. Section [TV] is devoted to the case of weak coupling, where very 
amazing effects are observed. It also contains numerical illustrations which played a crucial 
in our preliminary investigations. 

Next, in Sec. |Vj following the analogy with 3D, we question the validity of an "effective 
range" formula, and find sufficient (and sometimes necessary) conditions for the existence 
of the effective range. These conditions require a sufficiently-rapid average decrease of the 
potential at large distance r. We also find what happens when the potential decreases less 



rapidly, behaving like negative powers of r at large r. Finally, in Sec. |VI[ we return to 
the effective range in 3D, for the existence of which we believe to have the best sufficient 
conditions, in spite of the fact that the notion of effective range is known since a very long 
time and is well documented [71 El El QUI ITT] . 

II. THE CASE OF A FINITE RANGE POTENTIAL 

As an example, we discuss the case of a finite range potential 

V(r) = 0, for r>R, (2.1) 
The wave function for r > R, is given by 

*(Jfc,r) = C [J (kr) cos 8- Y (kr) sin 8] , (2.2) 
which has the asymptotic form 



* — ^ °o cos ( kr + S - • ( 2 - 3 ) 



We need to match at r — R the external wave function ( 2.2[ ) to the internal one, ^(k,r 



d r V(k, R) _k \J' {kR) sin 8 - Y^(kR) cos 8} 
V(k,R) ~ J {kR) sin 8 -Y {kR) cos 8 ' ^ ' ' 

where d r indicates derivation with respect to the second variable. 

By an extension of Poincare's theorem under conditions on V to be specified later, we 
have, if V(0,R) ^0 

d r ^(0,R) d r V(k,R) „ n2 



V(0,R) V{k,R) 



0(k 2 ) . (2.5) 
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On the other hand, for small z, the Bessel functions are given by 



J (z) = l + O(z 2 ) 



Y {z) = -Jo(z) 



71 



+ 0(z 2 



(2.6) 



where 7 is Euler's constant. 

The zero-energy wave function for r > R is given by 

\[>(0,r) oc In (r/a) . 



Combining Eqs.(2.4)-(2.7), we obtain 

2 



cot 5 



7T 



, 1 ka 



as k — > 



(2.7) 



(2.8) 



One should note that the 2D scattering length a gives the value of r where the asymptotic 
zero-energy wave function vanishes, in complete analogy to the 3D case. A major difference 
between the 2D and 3D cases is that, in 2D, the scattering length a is always non-negative. 



III. EXISTENCE OF THE SCATTERING LENGTH IN 2D 

In this section, we start with the S-wave (m = 0) solutions, u(k,r), which satisfy 

^ + ^ + k 2 - gV{r)\ u(k, r) = . (3.1) 



Of the two independent solutions of (3.1), under conditions on V(r) to be specified below, 
we choose as a regular solution 



u(k, 0) = , u(k, r) ~ \fr as r — > , 

corresponding to a wave function \l/(fc,r) = u(k,r)/ 'y/r which is finite at the origin. 
We introduce the Green's function G(r,r') for r, r' > 0, defined by 



(3.2) 



— — f- k 2 



G(r,' r) = 5(r — r') 



This G is given by 

G{r, r') 



71 



W [J (kr)Y (kr') - J (kr')Y (kr)] 9(r'-r) . 



(3.3) 



(3.4) 



The Volterra integral equation that gives the solution of Eq. (3.1) guaranteeing Eq. (3.2) is 



u(k,r) = u (k,r) + g / dr' G(r, r') V(r') u(k, r') 

'0 



(3.5) 



where we set u (k,r) = y/r J Q {kr 
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We proceed to find a bound on \G(r,r')\ which is uniform in k. We start with the 
well-known result 

Y (z) = -J (z) [ln(z/2) + 7] + R(z) , (3.6) 



7T 



where R(z) — > as z — > [12J. In Appendix [Aj we prove that for any z 



\R(z)\ < 



3tt 



Substituting Eqs. (3.6) and (3.7) into (3.4), we obtain 



\G(r,r')\ < \fr~r~ 1 In (J^j +3 



Now it is easy to show that 



where 



Hence 



ln + r 



In r if r > 1 
if r < 1 



|G(r,r')| < v^V ( v / 3 + 



In r 



In r' 



if r > 1 
• In r if r < 1 



In r' 



Using (3.11), we can deduce from (3.5) an integral inequality 



\u(r)\ < \fr + \fr ( V3 + — — I I(r) 



with 



I(r)=g J (V3 



In r' 



\V(r')\ \u(r')\ dr' . 



This inequality can be easily integrated, giving 

I(r) < exp 

and 



rr 

g r' (3 + hTr') \ V(r')\dr' 
Jo 



\u(r)\ < 2y/r (l + ln + r) cxp 
So, if the integral 

POO 

r (3 + I lnr|) |^(r)| dr 



g / r' (3 + I lnr'l) |V(r')| dr' 
. Jo 



converges, the iterative solution of (3.5) exists (see, for instance, [13 



(3.7) 



{3.1 



< In ^-J + 3 < ln + r + In" r + 3 < ( v 7 ^ + ( + — ^- J , (3.9) 



(3.10) 
(3.11) 

(3.12) 
(3.13) 

(3.14) 

(3.15) 
(3.16) 



We now proceed to derive an expression for the scattering length in 2D. Recalling the 
remarks in the introduction, the scattering length in 3D is given by the value of r for which 
the asymptotic reduced wave-function, i.e., ?i asym = A + Br, vanishes. In the 2D case, 



u 



asym 



(A + B lnr) y/r 



(3.17) 
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and we thus define the scattering length as 

a = exp [-A/B] 



(3.18) 



This shows explicitly that the scattering length is non negative. It is zero or infinity when 
5 = 0. 

The reduced radial equation in a partial wave reads — u" + (a/r 2 + V(r))u(r) = k 2 u(r) 
with a > —1/4. In the usual cases, a > —1/4, and the sign of the phase shift reflects 
how the effective potential departs from the pure centrifugal term. Here, a = —1/4 is the 
minimal admissible value. This is why, the departure due to V(r) seen in the low-energy 
phase-shift has always the same sign. 



The zero-energy limit of the integral equation (3.5), with Uo{r) = u(0,r), is 



u (r 



r + g 



/ rV In ^— J V(r') u (r') dr' 



(3.19) 



It can be rewritten as 
u (r) 



r 



9 



r' In ^— J V(r') uq{t') dr' — g I \fr' In ^— j V(r') Uo{r') dr' . (3.20) 



o 



We define X\ and X 2 as 



y/rV(r) uo(r) dr , 
r lnr V(r) Uo(r) dr 



(3.21) 
(3.22) 



Because of (3.15), valid for any k, these integrals converge. For the time being, we shall 
assume that X\(g) differs from zero. If X\(g) = 0, this constitutes the exceptional case 
discussed in Ref. PQ. 

From (3.20), we now get 



l-gX 2 (g)+gX 1 {g) lnr + g 



r' In — V(r') Uo(r') dr' 



Hence, as r — » oo, if we set 



In a 



9X 2 (g)-l 



(3.23) 



(3.24) 



gx 1 ( g ) ' 

the asymptotic uo(r) will have a zero at r = a. For any \g\ > 0, given our bound on uo(r), 
In a exists and is finite, and given by 



a = exp 



gx 2 



gXi 



(3.25) 



The next question is to see how this length a appears in the low-energy formula of 5(h). 
The asymptotic form of u(k,r) is 



lim u(k,r) = u asy (k,r) oc \Zr[J (kr) cos 5(k) —Yo(ky) sm6(k)] 



k fixed 



(3.26) 
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Hence from the integral equation for u(k,r) 



TT 



cot 5(k) = lim 



\--g\ Y (kr')V{r')u(k,r')Vr'dr' 



--g I J (kr')V(r')u(k,r')\/r'dr' 



(3.27) 



and from (3.6) we get 



1-9 



cot S(k) = - (ln(fc/2) + 7) + lim 

TT r^oo 



TT 



lnr' + — R(kr) V(r) u(k,r) y/r dr 



7f 

— —(7 / Jo(kr) V(r) u(k, r) y/r dr 



(3.28) 



Both integrals in (3.28) exist given our bound on \u(k,r)\ in Eq. (3.15) and the fact that 
R{kr')\ < 1, and if we also impose the second condition on the potential, namely 



ln 2 r |V(r)| r dr < 00, 



(3.29) 



we have 



cot 5(k) 



7T 



In 



+ 7 



1-„J 


/■oo 



TT 

lnr' + —R{kr) 


V(r) r) -y/r dr 




7T 

2 g J 


poo 

' Jo(kr)V(r)u(k,r) y/r dr 




(3.30) 



Next we show that \imk^ [cot 5(k) — (2/7r) ln(/c/2) + 7] exists. In (3.30), both integrals are 
of the type J °° F(k,r)dr where lim^Q F(k,r) exists for fixed r and \F(k,r)\ < B{r) with 
J °° B(r) dr < 00. Under this condition 



poo poo 

lim / F(k,r)dr= / F(0,r)dr. 
fc ^° Jo Jo 



(3.31) 



Hence we get 



fe^O 



lim < cot 5(k) 



In 



+ 7 



2 -1 + gX 2 _ 2 

7T fifXi 7T 



In a , 



(3.32) 



where X\ and X2 are given by Eqs. (3.21) and (3.22), and lna defined in Eq. (3.24). 
Thus we finally have 



k^O 



lim < cot 5(k) 



TT 



. 



(3.33) 



As shown before, this limit is uniform for small k. Hence as k — > 0, k > 0, we have 

2 



cot 5(k) 



TT 



hi 



ka 
Y 



+ 7 



o(l). 



(3.34) 
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This immediately leads to 



tan<5(/c) — 



tt/2 



oil) 



ln(fca/2)+7 |lnjfe| : 



(3.35) 



Thus under our general conditions on the potential V(r), the length a completely determines 
the (9(| \nk\~ 2 ) contribution to 5(k). 



Notice that our conditions (3.16) and (3.29), which can be summarized as 



POO 

/ (l + I In r I + (ln + r) 2 ) \V(r)\ r dr < 
Jo 



•x- 



(3.36) 



are in fact weaker than those used in [TJ. 

We now turn to th e case X\(g) = 0. Notice first that if X\ — 0, 1 — gX 2 ^ because 
otherwise, from (3.20) we would get uo{r)/^ — ► for r — > oo, which implies Uo(r) = 0. So 
for X\ = 0, uq oc yr, which corresponds to a zero-energy bound state and 



cot S(k) (ln(fc/2) + 7) -> 00 for k -> 



7T 



(3.37) 



If (? c is a value such that Xi(^ c ) = 0, X x has a single zero at g = g» c , because the successive 
zeros of X\ correspond to an increasing number of nodes for the zero-energy wave-function 
(this would not be true for the non-central case!). So, since 1 — g c X2(g c ) 7^ 0, we have 



either 



or 



00 




00 



for g 

for g 

for g 

for g 



9c g < 9c , 

9c 9 > 9c , 

9c g < 9c , 

9c 9> 9c ■ 



(3.38) 



IV. THE WEAK COUPLING CASE 

In this section, we shall reconcile the apparent contradiction between the sign of tan 5(k) 
as given by first order perturbation theory for small g, real k > 0, and on the other hand the 
universal sign for tan<5(fc) given by our results which for small k are independent of V(r). 

We will also explore further the behavior of the scattering length, a, as g — > 0. This turns 
out to be dependent on whether g J °° r V(r) dr is positive or negative. 



From Eq. (3.27), we have 



-g I J (kr')V{r')u(k,r')Vr'dr' 



tan<S(fc) = ^ , (4.1) 



7T 



g / Y (kr') V{r') u{k, r') vV dr' - 1 



The power series in g for u(k,r) is absolutely convergent. Thus for any real k > 0, tan 5(k) 
is given by the ratio of two entire functions in g. The perturbation series for tan5(fc) has a 
radius of convergence determined by the smallest zero, \go\ > 0, of the denominator. 



8 



Thus for small \g\, we have 



tan5(k) 



7T 



9 



r[J (kr)} 2 V{r) dr + 0(g 2 



(4.2) 



since u(k,r) = y/r J G (kr) + 0(g). 



On the other hand, we have from Eq. (|3.35|) 

tan5(fc) — 



tt/2 



0(1) 



ln(fca/2) +7 |lnfc| 2 ' 



(4.3) 



The sign of tan5(fc) in Eq. (4.2) changes when one goes from a repulsive gV(r) to an 
attractive force. To reconcile this with Eq. (3.35) or (4.3), we need to study further the 
dependence of a on g near g = 0. 



From Eq. (3.24), we have 



In a 



roc 

g / \fr lnr V(r) w (r) dr — 1 
Jo 



g I sfr V(r) u (r) dr 



(4.4) 



where Uo(r) satisfies the Volterra equation (3.19), which to first order in g reads 



Uo(r) = \fr + g \fr \ r' In — V{r')uo{r') dr' + 0(g ) . 



(4.5) 



Substituting this result in (4.4), we get for a 



a = exp 



oo oo r 

—l/g + J r\\irV(r)dr J rV(r)dr J r' V(r') \n — dr' 



+ 



r V(r) dr 



r V(r) dr 



0(g) 



(4.6) 



Hence if / °° r V(r) dr > 0, a(g) — > as g — > 0. However, if J °° r V(r) dr < 0, then as g — ► 
from gf > 0, a — > oo. 

In the first case, / °° r V(r) dr > 0, there is no bound state, tan 5 < 0, the scattering 
length a is small, and for small k, 5(k) ~ (tt/2)/ In fca. In Fig. [TJ is drawn for the 
potential gexp(-r), and several values of g. It is clearly observed that, while for k ~ 2 — 4, 
5(&) is nearly proportional to g, for small k there is an universal behavior. 



In this figure, the phase-shift is computed by two methods. First, the radial equation (3.1 ) 



is solved numerally and u(k,r) is matched into the Bessel functions, as per Eq. (3.26). For 
cross-checking, the radial equation is transformed into a non-linear, first order differential 
equation in r for the phase function 5(k, r) of Calogero |14| . and then 5(k) = lim^oo 5(k,r). 

For the second case J °° r V(r) dr < 0, we have at least one bound state for g > 0, even if 
g is small. However, the scattering length a is large for small g, and the asymptotic result 
5(k) ~ 1/lnfca will only be reached for very small k, k <^ 1/a. The phase-shift behaves as 
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FIG. 1: 8(k) for V = gexp(— r), and g = 
1/2, 1, 2. Each phase-shift is computed by two 
methods, but the two curves cannot be distin- 
guished. 



0.5 



1 1.5 



25 
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i 1 r 
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FIG. 2: 5{k) for V = -gexp(-r), and g = 1/2, 1, 2 (left) and g = 1, 2, 3 and 6 (right). 



shown in Fig. |2j For large tan<5(A;) > 0, in agreement with Eq. (4.2). An illustration is 
given in Fig. j2Tfor the potential V = —gexp(-r) and several values of the coupling g. 

Finally, one should note that Figs. [T] and [2] illustrate the Levinson theorem, whose validity 
in 3D was established by one of us, A.M. (and many others!), many years ago [3J HS1 HE], 
and this particular proof can be extended to 2D. The convention adopted in Fig. [2] is that 
S(k) — > as k — > oo. Hence 6(0) = nit, according to the number n of bound states. 



V. THE EFFECTIVE RANGE IN 2D 



Since we have been able to introduce the scattering length in 2D in analogy to 3D, it is 
natural to ask oneself if one can also define and effective range in 2D. In fact this has already 
been done by a group of Dutch physicists [TTj . However, what we want to do here is to find 
the conditions under which the effective range exists and what happens if it is infinite. We 



shall see that the condition (3.36) is not sufficient and that, crudely speaking, we need a 



potential decreasing faster at infinity. 

The method we shall use is, at least at the beginning, a carbon copy of the method used 
by Blatt and Jackson [18], Blatt and Weisskopf [7] and other textbooks [Sj. 

From now on, we shall use a different normalization for the wave function. First we define 
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the free zero-energy wave function Vo(r) which coincides with Uq{t) for r — > oo as 



vo{r) = \fr In (-J , 



(5.1) 



and the free wave function v(k,r) which coincides with the exact wave function at infinity 
will be normalized in such a way that 



v(k, r) — > v(0, r) = Vo(r) 



(5.2) 



for any finite r, as k — > 0, or more generally for any r(k) such that kr(k) — > for k — > 0. 
One such normalization is 



7T 



u(r) = [ co ^ ^ ^o(kr) — Y (kr)] . 



(5.3) 



On can check that in the limit k — > 0, the coefficient of y^r lnr is correct, and that cot 5 
disappears because of Eq. (3.34). 

The standard procedure [7] is to combine the exact and free equations for k = and for 
k 7^ to get 

/>oo 

lim[t>t>Q — v'vo] = k 2 {uuq — d^) dr , (5.4) 

r->0 /r, 



and substituting (|5.1|) and (|5.3 ) , 

- In 



cot S 



71 



ka 

y 



7 



k 



(vvq — uuq) dr . 



(5.5) 



This expression is exact. The crucial step is to say 

2 



lim — 

fe^o A; 2 



cot 5 — 



7T 



hi 



+ 7 



(ujj - ug) dr 



(5.6) 



where we use the fact that since v approaches vq and since the normalization of u is fixed 
by that of v , u approaches uq. However, there is absolutely no guarantee that the integral 
in the R.H.S. of (5.6) converges. 

Since uo approaches vq, it is sufficient to study (uq — vo)vo- Now, we use again a Volterra 
equation, but contrary to what was done in Sec. Ill, we start from infinity, i.e., 



uq = v + g V r 



In 



r' V(r') Uo(r') dr' 



(5.7) 



Since we only want to have the behavior of uq — Vq for large r, it is sufficient to take the 

(5.8) 



first term of the perturbative expansion of (5.7), i.e., 

^r'V(r') fo(r') dr' . 



u ~ v ~ g V?" 



Ill 



The integral in (|5.6|) will converge if 

v (r) y/rdr 



In 



r' V(r') Vq(t') dr' < oo 



(5.9) 
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and substituting (5.1) we get in the end 

V(r) r 3 



In 



dr < oo . 



(5.10) 



This is a sufficient condition for the existence of the effective range, defined by (5.6). 
However, if we restrict ourselves to potentials with a definite sign for large r, this condition 
is also necessary. If (5.9) diverges, the effective range is infinite. Hence (5.9) can be regarded 
as the best possible condition. 

Now, can we say something more when (5.9) diverges? We have investigated the case 
where 

V(r)=gr~ u for r>R with 2 < v < 4 . (5.11) 



Notice that v < 2 is incompatible with (3.36). The strategy is to write the R.H.S. of (5.5) 

as 



[(u - v)u + (uq - v )v] dr 



and use not only (5.7) but also the Volterra equation for u 



u = v + g \fr I [Jo(kr)Yo(kr') — J (kr')Y Q (kr)} v r' V{r') u{r') dr' 

J r 

and, since we are only interested in large r, we keep only 

POO 

u — v ~ g y/r / (J (kr)Y (kr') — J (kr')Y (kr)) vV V(r') v(r') dr' 



(5.12) 



(5.13) 



(5.14) 



Now comes the tedious, but straightforward task of estimating (5.12) by inserting (5.8) 



and (|5.14|), and the known expressions for v and vq. Details are given in Appendix |B| The 

ka s 



final answer is 
2 



cot 5(k) 



71 



111 



+ 7 



2 r2[ln( y f _±_ (±Y~ 3 T(u-2)T(2-u/2) 
7T [[n[ka)l v-2\2) T\v/2)T{v/2-l) 



, (5.15) 



if V = g r v for r > R. 



VI. EFFECTIVE RANGE EXPANSION IN 3D 

Even though the notion of effective range in 3D is well known (see, e.g., [TJ [8j [18] and 
references there to other pioneering contributions), we want to come back on the subject, 
because we believe that, at least to the best of our knowledge, its validity has not been 
investigated systematically. Mott and Massey [U] explain that if the potential decreases at 
large r as r~ s , the effective range will exist only for s > 5 1 . We shall give a somewhat 



broader sufficient condition, which could be the equivalent of (5.10) in 3D. K. Chadan has 



1 We thank the referee for a clarification on this point. 
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given to us a series of references on cases where the effective range formula is not valid [19], 
namely for a potential behaving like gr~ u at large r, with 3 < v < 5, one has 

k cot5 + - ~ g C(u) k u ~ 5 , (6.1) 
a 



where C(z/) is known. Again, this is analogous to Eq. ( |5.15 ) in 2D. 
Now, the effective-range formula in 3D is 

k cot 5 + - ~ \ r ok 2 , (6.2) 
for k — > 0, where the effective range is given by 

poo 

r = 2 (v 2 -u 2 )dr, (6.3) 
Jo 

(the factor 1/2 is such that r is about the radius in the case of square-well potential, see, 
e.g., [H])- Here v$ is the free reduced wave function at zero energy, with a normalization 

vo = 1 - - , (6.4) 
a 

and Uq is the exact solution normalized in such way that uq — Vq — » as r — > oo. So m 
satisfies the Volterra equation 

POO 

u (r) = v (r) + g / (r' — r) V(r')tto(r') dr' . (6.5) 



The asymptotic behavior of u — vq is therefore given by the lowest order iteration of (6.5) 

poo 

u — v o ~ / (r' — r)V(r')(l — r'/a)dr' , (6.6) 

J r 

and the convergence of (6.3) will depend on the convergence of the integral 

) POO 

(1 — r/a)dr / (r — r) V(r') (1 — r'/a) dr' . (6.7) 
The conclusion is that the effective-range formula will hold if 

POO 

/ r 4 \V(r)\ dr < oo . (6.8) 
Jo 

In addition, to get the existence of a scattering length, we should impose 

r | V(r) | dr < oo , 

POO 

/ r 2 | V(r) | dr < oo . 
Jo 



(6.9) 



If V has constant sign for large r, the condition (6.8) is not only sufficient, but also 



necessary. Then, if (6.8) diverges, 

k cot S(k) + l/a 



oo , (6.10) 



k 2 

for k — > 0. In particular, the existence of an effective range cannot be summarized as whether 
or not r 5 V(r) — > at large r. For instance 
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• For the potential 

such that r 5 V(r) —* 0, the effective-range formula is valid only if a > 1. 

• For the potential 

V(r) = exp [-r 12 sin 2 r] , (6.12) 

such that V{r) = 1 for r = nir, n£N, the effective range formula is valid, even though 
V(r) does not decreases faster than r~ 5 . In this case, it is better to think of (6.8) as 
a Lebesgue integral. 

A final remark: is generally believed to be positive, because at short distances, vq is 
much larger than Uq. In fact, this is not necessarily true. Take the potential 

It is easy to see that it has a scattering length a = 3R for any D > 0, and that the effective 
range goes to — oo for D —>■ oo, while it is positive for D close to 0. 



VII. OUTLOOK 

In this paper, we have presented a rigorous definition of the scattering length in two 
dimensions, both from the zero of the zero-energy solution and from the low-energy be- 
havior of the S'-wave phase-shift, and the equivalence between the two definitions has been 
demonstrated. The effective-range expansion is also given for short-range potentials, and 
conditions have been written for the existence of a scattering length and of an effective 
range. The effective-range expansion is also generalized for classes of power-law potentials 
with a longer range. 

The main result is that for short-range potentials, the S-wave phase-shift S(k) behaves 
at low energy such that 

2 1 f°° 

cotS=-[Hka/2)+ 1 } + -r k 2 + --- , r = 2 / (v 2 - u 2 ) dr , (7.1) 

where the zero-energy reduced wave- function Uq(t) has an asymtotic form Uo(r) — > vo(r), 
with Vo(r) being the free, zero-energy solution vq(t) = y/r ln(r/a). 

This investigation gave us the opportunity to revisit the same questions in three di- 
mensions, and to clarify and expand the existing results. An alternative derivation of our 
conditions (displayed in the first preprint version of this article) have been proposed recently 
by Chadan |H]. 

We believe that it is very unlikely that the existence of the scattering length in two or three 
dimension depends on the rotational symmetry of the potential V. In other words, under 
very general conditions on the potential V(r), without rotational symmetry, the scattering 
can be expected to be well defined. It is our intend to study this problem. 
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APPENDIX A: BOUND ON R(z) 



The function R(z) defined by Eqs. (3.6) has the following properties. 
1. From P2], 



R(z) 



71^ 



tt/2 



2. 



cos(z cos 0) ln(2 sin 0) d9 , 



/■tt/2 

77.(0) = , because / ln(2 sin 9) d9 = . 

Jo 



Proof: 



tt/2 i-k/2 i /.tt/2 

In(2sin0)d0= / ln(2 cos0) d0 = - / In(2sin20)d0 
o Jo 2 J 

In(2sin0)d0 = - / In (2 sin(0) d0 . 



o * 
which would lead to a contradiction unless it vanishes. 



o 



3. 



\R(z)\ < 



3tt 



(Al) 



(A2) 



(A3) 



(A4) 
(A5) 



Proof: 

8 r /2 i6 r /6 

\R(z)\ < — / |ln(2sin0)| d9 = — In 2 sin d0 , 
Jo 7T Jo 

But, due to the convexity of sin0 for e [0,7r/6], 2sin6 l > 9/(ir/6) in this interval, 
and the above integral is less than 7r/6, and hence 



|77(z)| < 



3tt 



(A6) 



On the other hand, introducing 



7(0) = / In(2sin0)d0 , 7(0) = /(tt/2) = , 7(0) < if e]0, 7r/2[ , (A7) 
o 

we get by integrating by parts 



R(z) 



Ti 



/•tt/2 

[cos(z cos 0) 7(0)]o /2 -z sin(z cos 0) sin 7(0) d0 

Jo 



(A8) 
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so 



\R(z)\ < 



2 fn/2 



7T- 



but by integrating again by parts 

-.71-/2 



sin# cos6I{6) d9 



/ sin 9 cos 9 1(9) d9 = - - / cos 29 ln(2 sin 9) d9 
Jo 4 J 

\R(z)\ _ 1 



so 



Also 



\R(z)\<^-, and Jim 2 



i r 

-Jo(z) + l = -l \1 -cos(z cos 9)} d9 



7T 



o z 2 2ix 

2 r 



(A9) 

(A10) 
(All) 



7T 



sin 2 [(zcosfl)/2] dfl 



< ^ r[(zcos9)/2] 2 d9= j 



(A12) 



APPENDIX B: DETAILS ABOUT THE EFFECTIVE RANGE IN 2D 

We want to estimate (5.12) in the case where J °°(vq — Mq) dr diverges. Since (15.12) is 
dominated by the large r behavior, it is sufficient to use the asymptotic form (|5.8|) or un — v n 



and (5.14) fo u — v. So we have 

POO 

/ [{v (u-v) + v (uo - v )] -X + Y 
Jo 

where 



(Bl) 



X = g- [ r In (-) dr f [J (kr)Y (kr') - J (kr')Y (kr)} r' 
2 J R \aJ J r 

x V{r') [cot S J (kr') - Y {kr')} dr' , (B2) 

and 

f°° f°° f r '\ i i ( r '\ 

Y = gj r [cot 5 Jo(kr) — Yo(kr)] dr J \n y— J V(r') r' \n y— j dr' 

or, by exchanging the order of integration, 



(B3) 



X = g\ 



71 



V(r) [cot 5 J (kr) - Y (kr)] dr 



and 



x J r' In H-J [J (kr)Y (kr') - J (kr')Y (kr)] dr' , (B4) 
r In Q V(r) dr J r In H-j [cot 5 J (kr') - Y (kr')] dr' . (B5) 
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If we restrict ourselves to V(r) = r v for r > R, we can use scaling, taking the variable 
z = kr, and get 

POO 

X = g-k»- 4 dzz 1 -" [cot 5 J (z)-Y (z)} 

Z JkR 

x jf z' In [ J (*) Y (z') - J (z')Y (z)} dz' , (B6) 



and 



Y = gk v ~ 4 [°° z x ~ v In ( dz f z' In f -) [cot 5 J (z') - Y (z')} dz' . (B7) 
Ai? ^kaJ J kR \z J 

Unless we have convergence problems at the origin, we can replace kR by zero for k — > 0. 
Retaining only the dominant terms, we get 



cot 5 



7T 



In 



ka 



+ 7 



-k 



v-2 



In 



ka 

Y 



1 2 



dz z 



7T 



dzV In - Jo CO 



+ 



dz^Jo(z) / d^^[J (/)W- ^o(^W)] 

Jo 



Now we have to evaluate certain integrals: 
1. From the Bessel differential equation 



2. 



3. 



/ z'J (z')dz' = -zJ' (z) , 
Jo 

f Z z%(z') dz' = -zY>(z) + - 

Jo 7T 

/ a: In f — J J (x)dx = 1 — Jo(-z) , 



Jo 

dz' = + M**o(*) = -^'(z) + - 

z->0 7T 

and using the Wronskian Yq(z)J (z) — Y (z)Jq(z) = 2/(irz), we get 
[ Z z'dz' [J (z')Y (z) - Jo(z)Y (z')} = - [1 - J (z)\ . 

JO 7T 

Overall, we get 



cot <5 



7T 



In 



Y 



+ 7 



i/-2 



~ g k 

= gk^ 2 



In 
In 



ka 
Y 

Y 



„i-f 



[1 - [Uz)f] dz 
2 2 f°° 

/ Jo^Ji^^dz 

n 2- is J 



(B8) 



(B9) 



(BIO) 



(Bll) 



(B12) 



(B13) 
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where the last integral can be evaluated as using a formula in [211 P- 35], [22, p. 715] and 
also [231 Vol. II, p. 52, Eq. (30)] where the conditions 2 < 3?ez/ < 4 are specified 

™™-~*-lM$" Zww\ (B14) 
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